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ABSTRACT 



The danger of unstable modes for space vehicles resulting from propellant 
sloshing is well known. Additional artificial stabilization through the control system can 
enhance the stability behavior of the vehicle if the gain values and sensor characteristics 
are properly chosen. Very useful results can be obtained from a simplified stability in- 
vestigation of an accelerometer controlled rigid space vehicle with respect to propellant 
sloshing. In this parametric study the liquid was simulated by an equivalent mechanical • 
model consisting either of raass-spring-damper systems or pendulum -damper systems 
for each sloshing mode and a nonsloshing mass with a moment of inertia. For simplifi- 
cation, only the propellant sloshing in one tank was considered. The influence of various 
parameters, such as ratio of sloshing propellant mass to total vehicle mass, natural 
slosh frequency, control frequency, control damping, gain factors of the attitude control 
system, accelerometer gain, its natural frequency, damping and location was investi- 
gated. The effect of elastic deformations, aerodynamic forces and the inertia of the 
swivel engines was not considered. Half of the thrust was available for control purposes. 
Only the dominant first sloshing mode was considered, since higher mode mass ratios 
are considerably smaller than the first sloshing mode (except in four quarter tanks, 
where the mass of one other sloshing mode is of considerable amount compared with that 
of the first sloshing mode) . 


The stability boundary was determined in terms of the amount of damping of the 
propellant in the tank required for various tank locations ( slosh mass location) . 

The potential hazard in control due to propellant sloshing in space vehicles can be 
eliminated by proper choice of tank form (slosh mass ratio decrease and slosh frequency 
increase) , proper selection of characteristics, location and gain value of the accelerom- 
eter, and as a last resort by baffles in the propellant tanks. ^ ( j 
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DEFINITION OF SYMBOLS 


m mass of spac e vehicle 

m sloshing mass 

s 

y lateral translation of rigid vehicle 

^s displacement of sloshing mass relative to the centerline of the tank 

F thrust 

angle of rotation of rigid vehicle relative to inertial space 

P gimbal angle of engines against centerline of the vehicle 

x coordinate of gimbal point (measured positive from center of gravity 

towards tail of the vehicle) 

g longitudinal acceleration of the vehicle 

x slosh mass location 

s 

I e ff moment of inertia of the vehicle about the center of gravity 

7 s slosh damping (twice the amount of damping factor) 

k 2 square of the radius of gyration 

k^ spring constant of slosh model 

o s natural circular sloshing frequency 

a tank radius 

h liquid height 

a 0 , a t gain values of attitude control system 

g 2 gain value of accelerometer control system 

indicated acceleration perpendicular to spacecraft axis 


IX 



DEFINITION OF SYMBOLS (Concl'd) 


coa 


u 

x a 

4 


* = gg2 


*CR 


accelerometer frequency ( radians/sec) 

accelerometer damping (amount of critical damping of accelerometer) 
accelerometer location 

ratio of sloshing mass to total vehicle mass 
gain parameter of accelerometer 
control frequency ( radians/sec) 
control damping 

location of center of instantaneous rotation 
radius of gyration without sloshing mass 


x 
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SUMMARY 


Due to the danger of unstable modes for space vehicles resulting from propellant 
sloshing, and since additional artificial stabilization through the control system can 
enhance the stability behavior of the vehicle, a simplified stability investigation was 
made of an accelerometer controlled rigid space vehicle with respect to propellant 
sloshing. 

In this parametric study, the liquid was simulated by an equivalent mechanical 
model consisting of mass-spring-damper systems, or by pendulum-damper systems 
for each sloshing mode and a nonsloshing mass with a moment of inertia. Only the 
propellant sloshing in one tank was considered, and the influence of various parameters 
was investigated. Half of the thrust was available for control purposes. Only the domi- 
nant first sloshing mode was considered. 

The stability boundary was determined in terms of the amount of damping of the 
propellant in the tank required for various tank locations. 

The potential hazard in control due to propellant sloshing can be eliminated by 
proper choice of tank form, proper selection of characteristics, location and gain value 
of the accelerometer, and as a last resor^by baffles in the propellant tanks. 


I. INTRODUCTION 


The dynamic response of the liquid propellant in the tanks of a space vehicle may 
lead to saturation of the control system. The alleviation of this influence can be achieved 
in many ways: proper choice of tank form, tank location, gain settings of the control 
system and, as a last resort, by the introduction of baffles. Additional artificial stabi- 
lization with angle -of- attack meter or accelerometer, which is used to alleviate the 



required control deflections of the gimbal engines, also has a strong influence on the 
stability of the vehicle due to sloshing. Very useful results can be obtained from a simplified 
stability investigation of an accelerometer controlled rigid space vehicle. The location 
and gain of an accelerometer and its vibrational characteristics, such as its natural fre- 
quency and damping, may enhance the stability character of the vehicle. In order to 
obtain a general knowledge of the influence of tank and accelerometer location, the inter- 
action between the translation, pitching, and propellant sloshing is investigated. The 
forces and moments of the oscillating propellant in the tank are described by an equivalent 
mechanical model, consisting of a mass-point attached to a spring and damper for each 
sloshing mode and a fixed mass with a moment of inertia (Fig. 1). Since the effect of 
higher modes is very small, only the first mode of sloshing is considered in the analysis. 
This is justified by the fact that, in a circular cylindrical tank, the ratio of the sloshing 
mass to the propellant mass for the second sloshing mode is already smaller than three 
per cent of the first sloshing mode. In a cylindrical tank with annular cross section, the 
second sloshing mode mass is, in the most unfavorable case of a diameter ratio of inner 
to outer diameter, less than 12 per cent. In the quarter tank, however, the mass of the 
second sloshing mode is about 43 per cent of that of the first mode. In this case, the 
second sloshing mode will have some influence and thus cannot be neglected. This case 
and the treatment of two-tank sloshing at various tank locations will be treated in another 
paper. 

For further simplification, forces and torques due to aerodynamic or elastic 
deformations of the space vehicle were neglected. The inertial effects of the swivel engines 
were not considered. The parameters affecting the stability are: the ratio of oscillating 
propellant mass to the total mass of the space vehicle, the natural slosh frequency, the 
control frequency, the control damping, the gain factors of the attitude control and accel- 
erometer control sensor and the location, the natural frequency and damping of the 
accelerometer. The stability boundary was determined in terms of the amount of damping 
of the propellant in the tank required for various tank locations (slosh mass location) . 

The effect of a real accelerometer was also investigated. In the cases where 
sloshing was suppressed, the stability of the system was represented by the gain value 
of the accelerometer versus its location. 

As an illustration of these effects on the stability, a vehicle was chosen whose 
dimensions are similar to the Saturn space vehicle. Four cases were treated: 

1. Stability analysis of an accelerometer controlled rigid spacecraft including 
the effect of sloshing and considering a real accelerometer with certain vibrational 
characteristics. 

2. Stability analysis of an accelerometer controlled rigid spacecraft including 
the effect of propellant sloshing and considering the accelerometer to be ideal. 
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3. Stability analysis of a rigid spacecraft with real accelerometer control (no 
sloshing included) . 

4. Stability analysis of a rigid spacecraft with ideal accelerometer control (no 
sloshing included) . 


II. EQUATIONS OF MOTION AND CONTROL EQUATION 


In the following investigation, all effects due to aerodynamics, structural defor- 
mations, and propellant flow through the tanks and pipelines (Coriolis effects) are neglect- 
ed. Furthermore, the inertial effects of the swiveling engines are not considered. A 
rigid space vehicle is treated, where the propellant in one tank is free to oscillate. The 
coordinate system has its origin in the center of gravity of the undisturbed space vehicle. 
The accelerated coordinate system is replaced by an inertial system so that the vehicle 
is* considered in an equivalent gravitational field. The x-coordinate of this inertial 
coordinate system is tangent to the standard flight path (Ref. 1). 

The curvature of the standard flight trajectory is neglected. Roll motion is 
neglected and the investigation is restricted to the interactions of translational and pitching 
motion and the propellant oscillations. For the sake of simplicity, the n disc motion/’ 
which describes the motion of the propellant in a completely filled and closed tank (solid- 
ified free propellant surface) , and the compliance of the gimbal engines are neglected. 
Furthermore, it is assumed that half of the thrust is available for control purposes. The 
equations of motion are then: 

1. Equation of Lateral Translation 

my + m s i s = F(cp+ ^P) (1) 

where m represents the total mass of the space vehicle, m g the sloshing mass of the first 
propellant mode and F the total thrust. The factor \ on the right hand side of Equation (1) 
is a result of only half of the thrust being available for gimbaling. The displacement 
oi the sloshing mass m s relative to the centerline oi the tank is represented by £s and the 
engine deflection by 0. 


2. Equation of Pitching Motion 
Fxf , 

I e ff 9 + 2 0 t m S x s - gm s £ s 

Ieff represents the effective moment of inertia of 
It can be written as: 


= 0 (2) 
the vehicle about its center of gravity. 
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Iq dx + m Q x§ + I 0 + tn s x| = mk 2 . 


(2a) 


^ef f 



dx + 


The distance of the swivel point of the gimbal engine from the origin is represented 
by X E> m A T tfr e mass per unit length of the air frame, I Q r its geometric moment of inertia, 
m 0 the fixed (nonsloshing) mass of the propellant, which is located at the station x 0 , Io is 
its moment of inertia about its center of gravity and Xg is the location at which the sloshing 
propellant mass is placed. The radius of gyration of the vehicle is represented by k. 


3. Equation of Propellant Motion 


m s!s + Hi 7 s 


m c 


6 s + k s £s " m s x s <P 


gm s cp + m s y = 0 


(3) 


t g 

uig is the natural circular frequency of the first propellant mode, “ is the damping factor 
of the propellant and k s is the spring constant of the equivalent mechanical model, describ- 

• _ ei . • i i • . ,i . i A i. 


ing the fluid motion in the tank /k 


"s _ . P 

— - Ujg 


section, the natural frequencies are 


For a cylindrical tank of circular cross 


‘n 


_ 

2 it 



8% 

a 


tanh 



(n = 1 , 2,3 ) 


(3a) 


where the radius of the tank is a, h is the propellant height and g is the longitudinal 
acceleration of the vehicle, en are the roots of the first derivative of the Bessel function 
of first order and first kind (Jj (en) =0) (Fig. 2). The mass m s of the sloshing 
propellant is obtained from: ( 2 , 3 ) 


“n = i"p 


2 tanh (eg |) 

(% l) - x ) 


(n = 1 , 2,3 ) 


(3b) 


where m p is the mass of the propellant in the tank (Fig. 4). 

For a cylindrical tank with annular cross section and diameter ratio k of inner-to- 
outer tank, the natural frequencies of the propellant are (Fig. 2a) 


1 /\ 

B s n - 1 

/ h \ 
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tanh 
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( 6 a- 1 I) 


(n = 1 , 2,3 ) 
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where ^n-i 


are the roots of the determinant ( Ref. 4) 


J' x a) (0 = q 

J' x (k£) Y i ( k ^) 

The mass ms of the sloshing propellant is obtained from ( Fig. 


n- 1 


m n- 1 


— m r 


2 

n §n-i 


k C x (k tanh j) 


(! - k b 


5 ) 


for n = 1, where m p is the mass of the propellant in the tank and 
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3t & 


n- l 


- k C x (k & 


O-l) 


* 2 l 1 


-V C l (k ^.3) (1 - k 2 §2.3) 


and 


c i ^n-a) 


^1 0^ 
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Yi (k Sn-i> 
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For four cylindrical quarter tanks with radius a and filled with liquid to a height h, the 
natural frequencies of the propellant are: (Fig. 3) 
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nm 


_1_ 

2 n 



%m 
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tanh 



(n = 1,2,3 ) 

(m = 0,1,2 ) 


where ^ are the roots of (^,) = 0. The mass of the sloshing propellant is 

obtained from 


64 - e nm tanh 


%n - 


(enm |) 


( c nm " ^ m2 ) Jfm ^ brnr^ (fnm a) U-° 


J 2m + 2 p. + 1 ( f -nm) 
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DO 


J 2m + 2p + 1 (e nm) 


(n = 1,2,3 ) 

(m = 0,1,2 ) 


M £2m( %m) 2 

(4m 2 - 1) € nm L, 

L (i=o 

for n = m = 1 where m p is the mass of the propellant in the tank. It can be seen that 
other modes cannot be neglected (Fig. 6) because, m 10 ~ 0.43 m n . (For m = 0 the fac- 
tor 64 is replaced by 32) . 

4. Control Equation 

P = a 0 cp + a i cp + g 2 

This equation represents the idealized control equation, relating the attitude angle q> 
and the indicated acceleration Ai of the accelerometer perpendicular to the space vehicle 
axis. Derivatives of (3 which produce increasing phase lags with increasing frequency 
have been neglected and are not important for the basic argument. The gain value of the 
accelerometer is represented by g 2 - 

5. Accelerometer Equation 

A* 2 £a : . 

-f + — — A i + Ai = y - x a cp - gcp 

<4a (5) 

where c% is the natural circular frequency of the accelerometer, is the damping 
of the accelerometer and x a is its location. If the accelerometer is considered ideal 
(w a » 1) its equation reads 

A i = y - X a Cp - gcp. (5a) 

By substitution of the control deflection p and the accelerometer equation into the 
equations of motion, the following set of linear homogeneous differential equations is 
obtained. 


g . . ~h - 

y - j (a 0 cp + a-L cp) - gcp + p i s - - A £ = 0 

<P + < a o 9 + a i cp> • ^ X S - ^2 g + 2p A i = 0 

y - x s cp - gcp + i 8 + ,7s Sa + ^ = 0 

+ = ^- A i + A i - y + * a $ + gcp = ° 

*4a 


( 6 ) 

(7) 

( 8 ) 
(9) 
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where 




m 


is the longitudinal acceleration 

is the ratio of the sloshing mass to the total vehicle mass 


is the distance of the swivel engine gimbal point to the center 
of gravity of the vehicle (origin) 

A = gg 2 is a gain parameter of the accelerometer. 


HI. STABILITY POLYNOMIAL 


Assuming solutions of the differential Equations 6, 7, 8, and 9 with the time 
dependency e , where s is a complex number, the stability polynominal can be obtained. 
The above differential equations are transformed into a set of linear homogeneous algebraic 
equations, of which the determinant of coefficients must be zero in order to have other 
than trivial solutions, i. e. 
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= 0 ( 10 ) 


i. Stability Analysis of Accelerometer Controlled Vehicle Including the Effects of 
Sloshing and Employing a Real Accelerometer 

Employing a real accelerometer with a natural circular frequency and a 
damping factor £ a results in the determinant (10). The evaluation of this determinant 
leads to a sixth degree polynomial. 

A 6 s 6 +A 5 s 5 +A 4 s 4 +A 3 s 3 +A 2 s 2 +A 1 s+A o =0 (11) 

of which the coefficients (A = 0,1,2, 6) can be written with 
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as the square of the circular control frequency (frequency of the pitch mode) and 
_ 2 a o 

a i " ut (13a > 

gX E a l = 2 r . , ( , _ A . Ax E x a \ 

2k 2 k ^ \ 2 2k 2 ) (13b) 


where £ c is the control damping and x a is the accelerometer location. 

A 6 = k 17 + k 16 X S 
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A 4 = k 9 + k 10 7 S + k i;L x s + k 12 Xg 
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A o = k o 
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(14b) 
(14c) 
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( 14e) 
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( 14g) 


where the abbreviations k A = 0,1, .. . 18) are functions of the slosh mass 

ratio, the control frequency, the control damping, the accelerometer gain value, its 
location and its vibrational characteristics, the slosh frequency, the attitude gain values, 
the distance from the gimbal point to the center of gravity of the vehicle and the radius 
of gyration of the space vehicle. Their respective values are: 
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The coefficients of the sixth degree polynomial (n = 6 ) are real. The nec- 
essary and sufficient conditions for stability are (Ref. 5): 

1. The coefficients A^, A^.-p ^.. 3 . .. >0 

{ A l3 A 0 > 0 if n is even! 

A o >0 if n is odd J 

2. The Hurwitz determinants 


H 3 > 0 if n is even 

Hn-i> *n -3 >0 J 

H 2 > 0 if n is odd 

This is for a sixth degree polynomial (n = 6 ) 

Ag> A 5 , A3, A 1? A 0 > 0 

h 5 , h 3 > 0 

which is 

A 5 a 3 A x 0 0 

Ag A 4 A ^ Aq 0 

H 5 = 0 A 5 A 3 A x 0 

0 a 6 a 4 a 2 a 0 

0 0 A3 A3 A 4 

A 5 A3 A-j_ 

H3 = Ag A 4 A 2 

0 A 5 A3 


(15) 

(16) 
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This means that all roots s = CT + i w have a negative real part. The conditions a = 0 
represent the dividing boundary between stability (damped oscillation) and instability 
(diverging oscillations). The boundary of the stability region is, therefore, 

H 5 = 0, A e = 0 (dynamic stability) (17) 

(A 0 = 0 static stability) 

It can be seen, that A 6 = 0 represents in the (x s , 7 S ) plane straight boundaries parallel 

to the 7 s - axis. It is x s = + — ^ k^. The Hurwitz determinant H 5 =0 can be repre- 
sented in form of ^ 

C 0 (x s ) + Cl (x s ) 7 S + c £ (x s ) 7|+ c 3 (x s ) 7s 3 + C 4 (x s ) 7 4 + c 5 ( Xs ) y 5 = o (18) 

and the (x s ) (A = 0,1, ... .5) are polynomials in xs. From Equation 18 the inter- 

section points with the x g - axis of the stability boundary curve can be obtained by taking 

Os = 0) and solving the equation for x : 

S 

c 0 (x s ) = 0 


Stability is achieved for all points (x g ,Y s ) above the boundary curve. The stability 
boundary breaks off at the left and right due to the results of A 6 = 0. Therefore, stability 
is only exhibited inside these boundaries, i. e. , above the stability curve, left of the right -- 
boundary of A@ = 0 and right of its left boundary. These latter boundaries are outside the 
space vehicle and are, therefore, due to the physical insignificance not indicated. 


2. Stability Analysis of an Accelerometer Controlled Vehicle Including the Effect 
of Sloshing and Employing an Ideal Accelerometer 

The employment of an ideal accelerometer (^1 » 1) instead of a real 


one with certain vibrational characteristics, 
a fourth degree polynomial 

A 4 s 4 + A3S 3 + Ags 2 + A x s + A 0 

of which the coefficients (A = Q, 1, 2, 
and 13 

A 0 = k 0 

A-i = kjL + k 2 7g 

A 2 = k 3 + k 4 7s + k 5 x s 
A 3 = k 6 + k y 7s + k 8 x s 
A 4 = k 9 + k 1Q Xg + k i;L xf 


results with Equations 1, 2, 3, 4, and 5a in 

: 0 (19) 

3, 4) can be written with Equations 12 

(20a) 

(20b) 

(20c) 

(20d) 

(20e) 
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where 
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The necessary and sufficient conditions for stability ( n = 4) are that 
1 • A 4? A 2> Aq ^ ^ 

2. H, >0 


which is 


= 


A 3 A 1 0 


*1 


- A 1 A 2 A 3 - AqA.§ - Af A 4 , 


The stability boundary is obtained from H 3 = 0 


A 4 = 0 


A jA ^A 3 — A oA 3 + A x 


( 21 ) 

( 22 ) 


which is (because the results are presented in the (x s , y s ) plane) 
(k x + k 2 x s + k 3 xf) + (k 4 +k 5 x s +k 6 xf) y s 
+ 7s (Ky + k 8 x s + k g Xg) + 7 s 0 


(23) 


with 


ki - k x k 3 k 6 - k x kg - k Q k 6 


k 2 k x k 5 kg + k x k 3 k 8 - 2k 0 k 6 k 8 - k x k xo 

_ _ 2 

k 3 = C-l k 5 k 8 - k Q k 8 - k x k xl 

k 4 = k x k 4 k 6 + k 2 k 3 k 6 + k x k 3 k 9 - 2k 0 k 6 k Y - 2k x k 2 k 9 

k R = ko L kc? + ki kcr k 7 + ki k 4 k R + k P k^ k fl - 2k n k 7 k R - 2k n k P k 


“2 t k 1 *5 *7 ' r *1 *4 *8 ^ ^2 *3 ^8 


O k 8 - tv 2 *io 


^6 ” ^2 ^5 ^8 “ 2k 1 k 2 k lx 
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ky = k 2 k 4 k^ ”k k k 4 ky "h k 2 k 3 ky ** ^0 ky — k 2 kg 

k q — k£ ky “k k 2 k 4 kg — k 2 k^Q 

k 9 = k 2 k 11 
k io = k 2 k 4 ky 


It can be seen that S 4 = 0 represents straight boundaries parallel to the 
7s - axis. From Equation 23 the intersection points of the stability boundary curve can 
be obtained by taking y s * 0 and solving the quadratic equation in x s . 

ki + k 2 x s + k^ Xg = 0 (24) 


Stability is achieved for all points (7 S > x s ) above the boundary curve. The 
stability boundary breaks off at the left and right due to the results of relation 21. There- 
fore, stability is only achieved inside these boundaries, i. e, , above the stability boundary 
due to Equation 23 and left of the right boundary of Equation 21 and right to the left bound- 
ary of Equation 21. Usually, the boundaries of Equation 21 are outside the space vehicle 
and are, therefore, not indicated, 

3. Stability Analysis with Real Accelerometer Control (Sloshing Not Included) 

The employment of a real accelerometer with the natural circular frequency 
and the damping factor results in the equations 


gcp - J p = 0 

(25) 

S X E 

— | p = 0 
2k 2 

(26) 

a o 9 + a i 9 + g 2 A i 

(27) 

2 U A i 

■ — + Aj_ = y - x a cp - gcp 

'^a 

(28) 


which are, with the introduction of the control Equation 27, 
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0 


( 29 ) 


y - gq) - | (a 0 cp + a x cp) - | Aj. = 

§ X E ✓ ■ 

9 + 2 p (a ° 9 + 31 ^ + ip A i = 0 ( 30 ) 

7 J + A i + A i - y + x a 9 + gcp = o ( 31 ) 

They can be obtained from Equations 6, 7 and 9 by setting the sloshing mass ratio M- = 0 
and leaving out the slosh Equation 8. 

Assuming again solutions of the form e st , where s = a + iw is a complex 
number, the stability polynomial can be obtained from the determinant. 


s 


2 


0 


-s 


2 


“ 2 ^ (2 + 3 .q + a ^ s) 
gx £ (a 0 + a x s) 

s*~ + 5 

2k 2 

x a s 2 + g 



if + Ik 

ui 


s + 


1 


= 0 


and is of the fourth degree 

S 4 s 4 + A 3 s 3 + A 2 s 2 + A.^ s + A q = 0 


(32) 


The stability boundary is presented such that the magnitude of gain of the acceler- 
ometer is determined versus the accelerometer location. Writing 

x a 

" a (33) 

X CR 


where |xcr | is the distance of the center of instantaneous rotation to the center of gravity 
of the space vehicle. The coefficients (A = 0, 1, 2, 3, 4) are then: 


Ao = k.Q + kj_A(l + 0) 

Ai = k 2 + k 3 A (1 + Q) 

A 2 = k 4 + k 5 A (1 + Q) 

A 3 = k 6 + ky A (1 + G) 



(34a) 

(34b) 

(34c) 

(34d) 

(34e) 
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where the abbreviations are the functions of the control frequency, the control damp- 
ing and the vibration characteristics of the accelerometer. 


k o ~ 


2 


k l = 


CJT 


k 2 ~ 


2 co c 


+ 2?c ut 


^7 ' “ 


£a 


oo = 




k 4 = 


X -}- ^ Ca £c 

oo a w a 


k, = - ^ - 


1 _ ^ 

2 2ojf ooa 


f _ ^ be (A b , 2 t a 
K 6 5 + 

oo a u, a 


k v - - 


£c °h 


<*Ja 


The stability boundary can finally be expressed in the Q) — plane by the 
equation 

ho + (1 + Q) + A 2 k 2 (l + Ct) 2 + A 3 k 3 (1 + Q) 3 = 0 ( 35 ) 

where 

ko = k ^ k 4 k 6 - k 0 k 6 - — - 

k i = k 3^4 k 6 + k 2 k 5 k 6 + k 2 k 4 k 7 ,- — 2 J * 3 - 2 k 0 k 6 t 7 - k^kf 

ooa 

k 2 ” ^3 k 5 k 6 k 3 k 4 k 7 k 2 k 5 k 7 “ — \ “ k O k 7 “ 2k ^ k^ ky 

tx) a 

k 3 - k 3 k 5 ky “ k-^ky 


It can be seen from k Q that the a - axis (A = 0) is in the stable region. 
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4. Stability Analysis with Ideal Accelerometer Control ( Sloshing Not Included) 


The case treating the stability with an ideal accelerometer can be obtained 


either by the Equations 25, 26, 27, and 5a or by working 
case. The stability polynomial is of second degree 


— = 0 
% 


with 


The stability boundary curve is obtained from A x = 0 and reads 

■k, 1 


A = 


k 3 1 + a 


in the previous 


s 2 + s + Aq — 0 

(36) 

= k 0 4- k x A (1 4* d) 

(37a) 

= k 2 + k 3 A (1 + Qj) 

(37b) 

= 1 - ^ (1 + 0) . 

(37c) 


and represents a hyperbola. 

The values k>, are 
k 0 = ujf 



k 2 = 2 £ c 
k 3 = -£c 

With these values the stability boundaries A x = A 2 = 0 are seen to be identical: 


A = . 

l + a 

The a - axis (A = 0)lies in the stable region. 


(38) 
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IV. APPLICATION 


The influence of the various physical parameters on the stability behavior of the 
vehicle is treated numerically by considering a Saturn type vehicle with the dimensions 
xe = 500 inches and a square of the radius of gyration of kg = 2. 513 x 10 5 (inch) 2 . The 
length of the vehicle was assumed to be 2300 inches with the center of gravity (origin of 
coordinate system) being 500 inches from the base of the vehicle. In those cases where 
the accelerometer location was not changed, it was chosen to be at x a = -240 inches, 
which is 240 inches forward of the center of gravity of the space vehicle. 

The parameters are: 

1. p * — ratio of sloshing mass to total mass of the vehicle. 

m 

2. Xg distance of swivel point from the center of gravity of the spacecraft. 

3. kg radius of gyration of the vehicle. 

4. u) undamped circular control frequency. 

5. control damping ratio. 

6. 10 s first natural undamped circular slosh frequency. 

7. ^s/2 damping factor of propellant. 

8. circular undamped frequency of the accelerometer. 

9. £ a accelerometer damping ratio. 

10. A = g g 2 product of longitudinal vehicle acceleration times accelerometer 

gain value g 2 . 

11. x & location of accelerometer measured from the center of gravity of 

the vehicle t a - (x cr is the center of 

J x cr| ' 

instantaneous rotation) . 

12. x c distance of the first sloshing mass from the center of gravity of the 

b 

vehicle. 

13. a Q gain value of the attitude control system. 
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Case A: 


The stability boundary of an accelerometer controlled vehicle with real or ideal 
accelerometer can be seen in this case. The motion of the propellant is suppressed. 
For an ideal accelerometer 1/u^ ~ 0 (u&» i) the stability is represented by the 
hyperbola 


1 + a 


where Oi 


is the ratio of the accelerometer location to the absolute value of 


'•cr I 


the center of instantaneous rotation. For real accelerometers, the stability boundary is 
not changing considerably, except to the amount that the stability region is slightly 
decreased; i. e. , the gain values A » g.g2 experience only stable values for smaller 
values than in the ideal accelerometer case. 


Table I. 


The magnitude of range of the parameters in Case A is given in the following 


TABLE I 
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Figure 7a 
Figure 7b 
Figure 7c 
Figure 7d 
Figure 7e 
Figure 7f 
Figure 7g 


It can be seen in Figure 7a that the ideal accelerometer is hardly different from 
the case of a real accelerometer with a circular natural frequency of 60 radians/sec. A 
decrease in the natural frequency of the accelerometer causes a slight shift of the stability 
boundary, thus decreasing the stable area slightly. Increase of the accelerometer damp- 
ing factor £ a from subcritical to high supercritical damping also decreases the stability 
area in the same sense as above (Figure 7b) for large accelerometer frequency (w a m 55 ), 

For small accelerometer frequency ( w a ■ 12.0) increasing accelerometer damping 
first increases the stability region slightly for very low damping values (until £ a ■ 0 . 2 ) 
then decreases the region again slightly with increasing damping £ a (Figure 7c) . In- 
creasing control frequency exhibits decreasing stability region for large and small 
accelerometer frequency (Figures 7d, e). For small accelerometer frequency, the 
decrease of the stability region is slightly more pronounced, since for small control 
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frequency ( in this case the stability boundary) it is already below the one with large 
accelerometer frequency. For changing control damping, the stability exhibits the same 
trends for both low and high accelerometer frequency. With increasing control damping, 
the stability region is slightly enlarged, especially noticeable in the case of small accel- 
erometer frequency. Thus the control damping has the strongest effect on the stability, 
especially for low accelerometer frequency. 

It may be noted, however, that due to other effects, such as bending, the accel- 
erometer should be located before the bending loop (where Y ? < 0) , In a rigid space 
vehicle the A value is also restricted to a certain range if aerodynamic forces are 
included in the treatment and the drift-minimum principle is considered. Without 
aerodynamic forces or for a neutrally stable rigid vehicle the drift-minimum principle 
is satisfied for A » a Q • This would mean in our case A * a Q * 3.5 . For 

smaller a Q values, a kind of load minimum (a 0 * 0) would be approached. (This effect 
of loads reduction, however, is only valid for quasi-steady states. ) Due to the drift- 
minimum consideration, the accelerometer is best located in the vicinity of the center 
of instantaneous rotation, thus having a larger range of values available. 

Case B: 


In the following Case B, the stability boundary with respect to propellant sloshing 
is treated for no accelerometer control. For ideal accelerometer ( w a » 1) , the 

stability equation was obtained in Section III, paragraph 2. If no accelerometer control 
is employed ( A ■ 0 ) , and only an attitude control system is used, the equations 

reduce to 

A 4 s 4 + A 3 s 3 + A 2 s 2 + Ai s + A Q * 0 

where the coefficients A v can be obtained from the equations 20a. . . . e and the k v 

from the by substituting for A * 0 . The boundaries due to A 4 * 0 are 



The stability boundary intersects the x s -axis in the values obtained from equation 
23 for 7 S = 0 and A = o and are for t^, 2 , 

2 ' 

o s 

x i = - | x cr I 


X2 = ^ 7 ^ 


(1 - H ) 


cr 
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For very small values of the stability boundary intersects the x Q -axis 

a 0 co s 2 J s 

close to the center of gravity (origin) and close to the center of instantaneous rotation. 

The value of x 2 is more sensitive to the changes of cu s , cu c and a Q . 

The magnitude of range of the parameters in Case B is given by Table II. 

If no accelerometer control is employed ( \ - 0) , and only an attitude control 
system is employed, it can be seen that propellant sloshing has a considerable effect on the 
stability of a rigid space vehicle. Figure 8a exhibits the danger zone of a rigid space ve- 

a > c 2 

hide for small values of — . It is approximately between the center of instantan- 

a o w s^ 

eous rotation and the center of gravity. With increasing slosh mass ratio p, the stability 
is decreasing in this region, i. e. , more baffling has to be employed in the tank in order 
to maintain stability. Considering only the wall friction of an unbaffled tank as contrib- 
uting to the damping (y =0. 02) , the vehicle is always stable for a slosh mass ratio of 
p < 0. 15. In Figure 8b the influence of an increase in the circular control frequency 
indicates a rather strong increase of the danger zone toward the tail of the craft. The 
sloshing mass ratio was considered to be p = 0. 1. The increase in the control frequency 
demands more baffling in the danger zone. For a control frequency of about twice the 
magnitude of the nominal control frequency of 0. 3 cycles/sec, about three times as much 
damping has to be introduced in the tank in order to maintain stability. The case of 
co c > a> s is very unfavorable since it demands even more baffling and exhibits a larger 
danger zone. Keeping the control frequency a considerable amount below the first 
sloshing mode frequency with relatively low damping in the tank avoids instability, and 
keeps the danger zone restricted to the approximate region between the center of gravity 
and the center of instantaneous rotation. For a sloshing mass ratio of 10% the wall 
friction is sufficient to maintain stability as long as the control frequency is below 0. 3 
cycles. The change of the control damping can be seen in Figure 8c. It exhibits 
the following trend: for increasing subcritical control damping (f c <1) , the stability 
is decreasing in the danger zone, i.e. , more damping has to be introduced to maintain 
stability. The danger zone is unchanged by the change of the control damping. For in- 
creasing supercritical damping (f c > 1) the stability is increasing, i.e. , less damping is 
required in tank locations in the danger zone. No baffling is required for a sloshing mass 
of p = 0. 1 in the danger zone, if the control damping £ c < 0. 5 or if £ c < 2. 5. This means 
that for the parameters considered (cu s = 5.0, o> c = 2, a Q = 3,5, simple attitude control sys- 
tems with lead net work) ,the wall friction in the tank is sufficient to maintain stability for a 
rigid space vehicle of the Saturn type. A further question that always presents itself to the 
design engineer is the problem of tank design. Tanks with large diameters exhibit low sloshing 
frequencies which are in many cases too close to the control frequency of the vehicle. Provid- 
ing partitions in the tank or clustering the tanks increases not only the sloshing frequency 
but also reduces the sloshing mass considerably. In Figure 8d the effect of the change 
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of the sloshing frequency is shown. The danger zone increases toward the aft of the ve- 
hicle for decreasing slosh frequency which decreases the stability considerably. An 
increase of the slosh frequency enhances the stability, i. e. , decreases the danger zone 
and reduces the amount of baffling necessary to maintain stability. For the treated case 
of a rigid space vehicle with a sloshing mass ratio of P = 0.1 > the wall friction will be 
sufficient to maintain stability at any tank location if the natural sloshing frequency is 
above 3. 0 cycles/ sec. 

The control factor a Q of the attitude control systems exhibits a slight decrease of 
the danger zone toward the base of the craft and a slight decrease in stability for decreas- 
ing values (Figure 8e). 

Case C: 


The stability boundaries with respect to propellant sloshing are treated for an 
accelerometer controlled spacecraft employing an ideal accelerometer ( w a» 1) in 
addition to the control system of the previous Case B. From this it can be seen that the 
left and right boundaries depend strongly on the accelerometer locations x a and the 
accelerometer gain values A . The mass ratio is also of some importance, especially 
for the left boundary, i. e. , for the negative sign in front of the square root. For most of 
the treated cases, these boundaries are beyond the spacecraft and therefore not indicated. 

Figure 9a shows the influence of increasing slosh mass ratio. The most favorable 
gain value A * 1.0 for an ideal accelerometer controlled space vehicle was employed. 
The danger area is decreased considerably to a very short zone below the center of 
instantaneous rotation where wall friction for the liquid in the tanks already is sufficient 
to maintain stability. The increase of the control frequency above the slosh frequency 
results in a decrease in stability and an increase of the danger zone towards the base of 
the vehicle. Due to ideal accelerometer control , however, the damping necessary to obtain 
stability for sloshing mass of 10 of the vehicle mass is very small, even in the case when 
the control frequency is twice the sloshing frequency. For * 5.0 radians/sec and 
wc =* 10.0 radians/sec, the necessary damping As “ 0.01 (Figure 9b). The change 

of the control damping exhibits the same behavior as in the previous case. It decreases 
the stability for increasing subcritical damping and increases the stability for increasing 
supercritical damping. The effect in the already small danger zone is hardly recognized, 
(Figure 9c) . If the slosh frequency is below the control frequency w s < w c , the 
danger zone is large and covers nearly the complete region behind the center of instan- 
taneous rotation. Small damping, however, already guarantees stability (Figure 9d) . 

For increasing slosh frequency, a considerable decrease in the danger zone is noticed 
and an increase in stability is obtained. Figure 9c shows that no appreciable stability 
is encountered by change of the attitude control value a Q . The accelerometer gain value 
A = g 2 g has a very pronounced influence on the stability behavior of the vehicle. In 
Figure 9f it can be seen that, for A < 1 , the danger zone is behind the center 

of instantaneous rotation and increases toward a value slightly behind the center of gravity 
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as A approaches the value z.ero. A loss in stability takes place in the danger zone 
while A is decreasing. For increasing A > 1, the danger zone shifts forward of 
the center of instantaneous rotation and increases the magnitude of the zone and the 
required damping for tanks in this location. From a certain value A > 1.5 on, 
propellant motion in ak *^nks in front of the center of instantaneous rotation leads to 
instabilities while, for tanks behind the center of instantaneous rotation, the motion 
of the propellant is stable. The value A - l is the most favorable gain for which 
all other parameter changes were performed. In Figure 9g, the influence of the location 
of the accelerometer upon the stability can be seen. For the most favorable gain value 
A ■ 1, th e change of the location of the accelerometer did not lead to instabilities. For 
other gain values, however, the location of the accelerometer definitely has an influence 
on the stability. For gain values A larger than 1. 5 , the fluid motion in all tanks behind 
the center of instantaneous rotation exhibits strong instabilities for accelerometer loca- 
tions behind the center of gravity. The propellant motion in' tanks in front of the center 
of instantaneous rotation shows instabilities if the accelerometer is located in front of 
the center of gravity (Figure 9h) , For gain values A smaller than unity, the behavior 
of the stability boundary is similar to that one without accelerometer control. The pre-r 
vious results are only valid for an ideal accelerometer, which means for an accelerometer 
frequency of large value ( ^ » 1 ) . For an accelerometer with vibrational character- 
istics w a and £a (noted as a real accelerometer) the stability behavior is quite different 

from that of an ideal accelerometer depending mainly on the value of its natural frequency. 


The range of parameter changes can be depicted from Table II. Two circular 
frequencies for the accelerometer were treated: co a = 12 and 55 radians/sec. In Figure 
10a, the influence of increasing sloshing mass ratio can be seen for a space vehicle with 
additional accelerometer control, where the accelerometer has a circular natural fre- 
quency of cc a = 55 radians/sec and a subcritical damping of £ a = 0. 7. With increasing 
slosh mass ratio, the danger zone increases aft of the center of instantaneous rotation 
and the stability decreases. That is, in the danger zone more damping is needed to main- 
tain stability. It can be noted, however, that for wall friction only (y s = 0. 02) , the ve- 
hicle is stable for sloshing mass ratio p < 0. 30. Increasing the gain value to X = 1. 5 
(an optimal value for stability) with an accelerometer of circular natural frequency of 
ca a = 55 radians/sec, no instabilities are obtained for slosh mass ratios of p< 0. 21 (Fig- 
ure 10b). For larger slosh mass ratios the vehicle becomes very strongly instable. For 
low accelerometer frequency (u; a =? 12 radians/sec) , the stability situation becomes 
hopeless from the standpoint of sloshing ( Figure 10c) . Sloshing is even excited by the 
accelerometer control and creates a situation which is worse than in a vehicle with simple 
attitude control only. The danger zone is increased and continues to increase with in- 
creasing slosh mass ratio. Very high bafflingwould be necessary in orderto maintain sta- 
bility. The damping required in this case is about three to four times larger than in the case 
without accelerometer control. From this we can conclude that the accelerometer fre- 
quency should be well apart from the slosh control frequency. With increasing gain value 
A. this situation becomes even more unfavorable; i.e. , the danger zone increases over 
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almost the entire vehicle and stability requires tremendous damping values. Figure lOd 
exhibits the effect of the change of the control frequency for a vehicle with an accelerom- 
eter frequency of w a “ 55 radians/ sec. An increase of the control frequency below the 
slosh frequency decreases me danger zone slightly, while an increase of the control fre- 
quency above the slosh frequency increases the danger zone towards the aft of the vehicle. 
The damping during the increase of the control frequency remains comparatively small 
( 7 a “ 0.01 and less) . For an accelerometer frequency of w a - 12 radians/sec, 

the stability behavior can be seen in Figure 10c. Increase of the control frequency 
results in a decrease of the danger zone and an increase in stability. The effect is more 
pronounced than in the previous case, and the vehicle is less stable for equal control fre- 
quencies in comparison with the case of w a . 55 radians/sec. Figure lOf shows the 
influence of control damping. Increasing control damping decreases the danger zone and 
increases the stability. For the use of the low accelerometer frequency Hi >■ 12 
radians/sec, the stability behavior can be seen in Figure lOg, which exhibits, for increas- 
ing subcritical damping, decreasing danger zone and increasing stability, while for 
increasing supercritical damping increase of the danger zone and loss of stability area 
are encountered. Low subcritical damping, however, seems unfavorable in both acceler- 
ometer frequency cases; especially in the low accelerometer frequency case, more 
baffling has to be applied in order to maintain stability. The effect of changes in the 
sloshing frequency is shown in Figures iOh and iOi. Slosh frequencies below the control 
frequency exhibit a large danger zone from a value slightly in front of the center of instan- 
taneous rotation almost to the base of the space vehicle. For low accelerometer frequency 
(w a '» 12 radlans/sec) > this is even more pronounced. Increasing slosh frequency shows 
a decrease in the danger zone and an increase in stability. Increasing slosh frequency 
above the control frequency w c <Ct, s <0J a decreases the stability and increases the danger 
zone towards the base of the vehicle again. For low accelerometer frequencies, these effects 
are more magnified. For a slosh frequency in the vicinity of the accelerometer frequency, 
the stability is decreased and the danger zone is increased. The gain value a Q of the atti- 
tude control system has only very small influence upon the stability boundary for an accel- 
erometer frequency co a m 55 radians/ sec which is large compared to the control 
frequency w c “ 2 radians/sec (Figure iOj). The danger zone is only a short 

length behind the center of instantaneous rotation. With the wall friction damping value 
( 7 a * 0. 02) , the vehicle is already stable for the given parameters. For small 
accelerometer frequency (w a = 12 radians/sec) , a small effect can be observed. The 
danger zone and stability region, however, are unfavorable (Figure 10k). 

Figures 10 1 and 10 m exhibit the influence of the change of the accelerometer 
gain value ^ “ 8s 8 • For an accelerometer frequency of w a = 55 radians/sec, an 

increase of the gain value to about A* 1.5 enhances the stability and decreases the danger 
zone to a very small region between center of instantaneous rotation and the center of 
gravity. Further increase of A extends the danger zone in front of the center of instan- 
taneous rotation. The larger A becomes, the more baffling is required. If the acceler- 
ometer frequency is onlyo) a - 12 radians/sec, the situation is quite different. For 
increasing gain value, the danger zone increases from the center of instantaneous rotation 
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to the aft of the vehicle. Very strong dam ping in the tanks has to be introduced in order to 
maintain stability. Further increases decrease the stability region again, but still require 
a large zone to be baffled very strongly , while in the case of a large accelerometer, low 
damping is already sufficient for stability. A damping twenty to thirty times larger is 
required for the case of low accelerometer frequency. Again we can conclude that large 
accelerometer frequency eliminates the greatest part of the problem. To obtain the 
influence of the very important parameter, the natural frequency of the accelerometer, 
Figure lOn, shows that for increasing accelerometer frequency, the danger zone enlarges 
and requires more damping. Above the slosh frequency for >2 > a decrease 

of the danger zone and an increase in the stability area are noticed. The larger the ratio 

cj a 

— the less damping is required to obtain stability. 

The increase of the damping of the accelerometer exhibits an increase of the 
danger zone and requires more damping in the tanks. This effect is more magnified for 
small accelerometer frequencies (Figures 10 o, 10 p, and 10 q) . It can be noticed that 
from twice the critical damping on, for small accelerometer frequency w a * 12 radians/ 
sec, a further increase of the accelerometer damping slightly decreases the danger zone 
and slightly enhances the stability. A very important parameter in the design of a control 
system is the location of an accelerometer as indicated in Figure lOr, 10s, and lot. A 
location of the accelerometer behind the center of gravity must be avoided. For a natural 
frequency of = 55 radians/sec of the accelerometer, any location in front of the 
center of gravity requires low damping only for stability in a small danger zone. For an 
acceleration frequency w a « 12 radians/sec, a location of an accelerometer behind the 
center of gravity requires high baffling in the tanks between the base and the center of 
instantaneous rotation. Further shifting of the accelerometer toward the nose of the 
vehicle enhances the stability and decreases the danger zone. All these results are valid 
for a rigid space vehicle. It should , however, be noted that bending in an accelerometer- 
controlled vehicle plays an important role. Gain factor A , location of the accelerometer, 
and the ratio of the bending frequency to the accelerometer frequency as well as to the 
control frequency have some definite influence. If the bending and control frequency are 
sufficiently apart, only negative bending displacements are permissible. This indicates 
that for the control of the first two bending modes a location in front of the center of 
gravity, where both bending modes exhibit negative deflection, is favorable. This location 
also would remedy the sloshing problems with little baffling, but since it exhibits limitations 
in magnitude due to rigid body control (Case A) , it must be treated with care with respect 
to gain value A of the accelerometer. 

V. CONCLUSION 


For space vehicles of increasing size, the influence of propellant sloshing upon the 
stability becomes more pronounced. Aerodynamically unstable space vehicles create loads 
and require control torques, which can be a potential hazard to the flight performance of 
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the vehicle. Therefore, artificial stabilization such as angle-of-attack meter and/or 
accelerometer control is employed. This helps to alleviate the required control deflec- 
tions of the gimbal engines, jet vanes, or air vanes. For a given configuration, the control 
requirements are strongly dependent upon the gain settings of the control system (attitude 
and accelerometer control). Sloshing propellant requires additional control torques. In 
order to minimize this effect, proper choice of the control values, tank forms, tank 
location, and, as a last resort, baffles will enhance the stability situation. The amount 
of sloshing mass is determined by the tank geometry, the mass density of the propellant, 
and the liquid height. The tank location, i. e. , slosh mass location, plays an important 
role. The conclusions drawn from the results of a rigid space vehicle with additional 
ideal accelerometer control lead to optimistic results concerning the baffling of the tanks 
in order to maintain stability. An actual accelerometer has its own natural frequency; 
therefore, the relation of control, sloshing and accelerometer frequency are of utmost 
importance. Without consideration of propellant sloshing, rigid body stability can be 
obtained for any location of the accelerometer along the vehicle. Drift minimum require- 
ments, however, demand certain values for the gain ^ “(g^g) > which can only be 
obtained in front of the instantaneous center of rotation. However, for bending mode 
stability, the accelerometer should be located between the front nodal point and the center 
of instantaneous rotation for bending frequency of at least three to five times the control 
’frequency. In most cases, a simple control system such as that considered here in the 
rigid body analysis cannot be employed in bending feedback analysis without additional 
phase shaping in the network. 

For w c < °°s << w a the amount of damping required in the tanks to main- 
tain stability decreases considerably. One further advantage of an accelerometer controlled 
vehicle is the decrease of the danger zone, i. e. , the zone where, for a sloshing mass 
being located, a baffle of appropriate damping characteristic must be introduced to obtain 

stability. This danger zone was for values < i for a rigid vehicle with 

gj 2 a o 

simple attitude control system with lead network between the center of gravity and the 
center of instantaneous rotation. With additional accelerometer control, this danger zone 
can be decreased if the accelerometer frequency is large enough (co a » oj c ) . For small 
accelerometer frequency (w a 5 to 10 •u c ) the danger zone is enlarged for 

decreasing ,-u a . 

In the following the effect of the various parameters will be summarized: 

1. Slosh Mass Ratio P 

a. Without Accelerometer Control : Increasing slosh mass ratio^ecreases the 

__ L °c ' 

stability and requires more baffling. The danger zone is for values “ " < 1 

d o 

approximately between the center of gravity and the center of instantaneous rotation 
(Figure 8a) . 
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b. With Accelerometer Control 


((1) w a - 55 radians/sec: Increasing slosh mass ratio decreases the 
stability and requires more baffling. The damping required to maintain stability is only 
two-fifths of that of a rigid vehicle without accelerometer control. The danger zone is 
decreased to three-fifths of the danger zone for a vehicle without accelerometer control 
(starting at the center of instantaneous rotation toward the aft of the vehicle. Figure 10a) . 

( (2) w a _ 12 radians/sec: Increasing slosh mass ratio decreases the 
stability and requires more baffling. Ten times as much damping is needed in the tanks 
compared with the case of accelerometer frequency ^ ■ 55 radians/ sec. Even five 
times as much damping has to be introduced in the tanks compared with the results of a 
space vehicle without accelerometer control. Thus it can be concluded that no acceler- 
ometer control should be employed rather than one with too low a natural accelerometer 
frequency. The danger zone is even larger than for a system without accelerometer 
control and increases toward the aft of the vehicle for increasing slosh mass ratio (Figure 
10c) . 


(3) Ideal accelerometer: An ideal accelerometer eliminates the danger 
zone almost completely except for a small region at the center of instantaneous rotation, 
where wall friction is already sufficient to maintain stability. 

2. Control Frequency w c 

a. Without Accelerometer Control : Increase of the control frequency decreases 
the stability and increases the danger zone toward the aft of the vehicle. Furthermore, it 
requires stronger baffling. The situation becomes worse for control frequencies larger 
than the slosh frequencies ( Figure 8b) . 

b. With Accelerometer Control 

(1) w a m 55 radians/sec: Increasing control frequency decreases the 
stability slightly and increases the danger zone toward the base of the vehicle. The baffle 
requirements are very low compared to the case without accelerometer. Figure lOd shows 
that wall friction is already sufficient to maintain stability. 

( 2) w a m ^ radians/ sec: Decrease of the control frequency decreases 
the stability and also increases the danger zone toward the aft of the vehicle (Figure lOe). 

( 3) Ideal accelerometer: Increase of the control frequency decreases 
the stability and increases the danger zone toward the base of the vehicle ( Figure 9b) . 

Again as in the cases 1 and 2 , the baffle requirements are smaller than in the case 
without accelerometer control. Except for a control frequency of a w a _ 12 radians/ 
sec, the same baffling is required in the case without accelerometer control lor 

w c ■ 10 radiftns/sec. 


29 


3. Control Damping £c 


a. Without Accelerometer Control : Increasing subcritical control damping 
decreases stability, whereas increasing supercritical damping increases the stability. 
Supercritical control damping requires less baffling ( Figure 8 c) . 

b. With Accelerometer Control 

< 1 ) U) a > 55 radians/sec: Increasing subcritical damping increases 
the stability and decreases the danger zone. Increasing supercritical damping increases 
the danger zone and decreases the stability slightly. The baffling required in the tanks 
is very small for all control damping values. It is of the same magnitude as in the case 
without accelerometer control (Figure lOf), 

( 2 ) = 12 radians/sec: Increasing subcritical control damping 
increases the stability and decreases the danger zone. An increase in supercritical 
damping first increases the stability and decreases the danger zone, but around the 
propellant sloshing frequency, it decreases the stability and increases the danger zone 
toward the aft of the vehicle. The damping values required for stability are at least two 
to three times larger than those of the two previous cases (Figure lOg). 

(3) Ideal accelerometer: For the particular gain value of A * 1 the 
change of control damping has practically no influence. The trend, however, is 
the same as in the case witn no accelerometer control (Figure 9c). 

4, Sloshing Frequency 

a. Without Accelerometer Control : Increasing slosh frequency enhances the 
stability and decreases the danger zone. For slosh frequencies below the control frequency 
the danger zone is the entire length behind the center of instantaneous rotation and requires 
considerably more damping to maintain stability (Figure 8 d) . 

b. With Accelerometer Control 

(1) w a =■ 55 radians/ sec: For the slosh frequency below the control 

frequency, the danger zone is enlarged and the stability is decreased. For increasing 
slosh frequency, the danger zone first decreases, and the stability increases slightly. 
Further increase in the sloshing frequency toward the accelerometer frequency decreases 
the stability and increases the danger zone (Figure lOh). 

(2) oo a = 12 radians/sec: This case exhibits the same behavior (Figure 
10 ) except for about three times larger damping being required to maintain stability. 
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(3) Ideal accelerometer: Increasing the sloshing frequency slightly 
decreases stability. Sloshing frequency below the control frequency increases the danger 
zone and decreases the stability (Figure 9d) . 

5. Gain Value a 0 of the Attitude Control System 


a. Without Accelerometer Control : A decrease in the gain value ao decreases 
the stability and increases the danger zone. It is, however, of slight effect (Figure 8e) . 

b. With Accelerometer Control 

(1) w a = 55 radians/sec: Change of a Q exhibits practically no effect 

(Figure lOj) . 

(2) UJ a = 12 radians/sec: The change of a 0 also shows only minor effects. 
Increasing a Q decreases the danger zone slightly and causes a minor increase in stability. 
The danger zone was already three times as large as in case 1 and twice as large as in 
the case without accelerometer control. The damping required is about ten times that 

of the case 1 and about twice as much as in the case with no accelerometer control. 

( 3) Ideal accelerometer: The change of a Q has practically no influence. 

6. Gain Value, ^ ~ 9 of Accelerometer 

a. >a ^ radians/sec: For increasing gain value A , the danger zone 
between the center of instantaneous rotation and the center of gravity decreases accom- 
panied by increasing stability. At A *=1.5 the optimum seems to be reached. For 
further increasing gain value , the danger zone is in front of the center of instantaneous 
rotation and requires increased baffling (Figure 10k). 

b. radians/ sec : For increasing gain value A the danger zone 

increases aft of the center of instantaneous rotation and requires excessively large 
baffling. Any accelerometer control with a low natural frequency has detrimental effects 
upon the stability (Figure 10m). 

c. Ideal accelerometer : For increasing gain value A , the danger zone 
between the center of instantaneous rotation and the center of gravity of the vehicle is 
decreased and the stability is increased. Here a = l.Q seems to be the optimum 
value, while for further increase of A , the danger zone increases in front of the center 
of instantaneous rotation and requires more damping to maintain stability (Figure 9f) . 
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7. Accelerometer Frequency w a 


Increasing accelerometer frequency decreases the stability and increases 

the danger zone behind the center of instantaneous rotation for values of -cj** < | , 

i. e. , for values of the accelerometer frequency below the sloshing frequency. For 
further increase of the accelerometer frequency, the stability enhances and the danger 
zone decreases. For ideal accelerometer control ( w a — > °° ) , the danger zone is 

diminished to a small region about the center of instantaneous rotation, where practically 
no baffling is required to maintain stability (Figure lOn). 


8. Accelerometer Damping £a_ 


a. radians/sec: An increase of accelerometer damping ^ 

results in an increase of the danger zone aft of the center of instantaneous rotation and 
a decrease in stability (Figure 10 o). 

b. ^ a ~ 12 radians/sec: Increase of the accelerometer damping increases 

the danger zone aft of the center of instantaneous rotation and decreases the stability for 
subcritical and low supercritical damping. For supercritical damping £ a > 2 , further 
increase in damping has the opposite effect. The danger zone for corresponding damping 
values is larger than in the previous case oo a ■ 55 radians/ sec for £ a < 3 and vice 
versa for £ a > 3 . For low subcritical damping, the required baffling is approx- 
imately the same in both cases (w a = 55 and u a = 12 radians/sec), w hile for large sub- 
critical damping, twice as much damping is needed in the tanks to maintain stability for 
the low accelerometer case. For low supercritical accelerometer damping, about four 
times as much baffling is needed, while for large supercritical damping only about twice 
as much baffling is required in the tanks compared to Case A (Figures lOp and lOq) . 

f x a 

9. Accelerometer Locations x a ora * 

V J. X CR I 

a. = 55 radians/sec: Accelerometer location behind the center of gravity 
has a detrimental effect on the stability and exhibits only a rather short region of stability 
immediately in front of the center of instantaneous rotation. This stability region increases 
toward the nose of the craft as the accelerometer location is moved from the aft of the 
vehicle toward the center of gravity. Accelerometer locations in front of the center of 
gravity exhibit only small danger zones between the center of instantaneous rotation and 
the center of gravity, where only slight baffling is required to maintain stability. (Figure 
lOn, wall friction is in many cases already sufficient. ) 
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b. = 12 radians/sec: Again, as in the low accelerometer frequency 
case, stability is endangered. Location of the accelerometer behind the center of gravity 
again exhibits detrimental effects upon the stability and results in a danger zone which 
reaches from the center of instantaneous rotation to the tail of the vehicle. Large baffling 
would be required to maintain stability. But even for an accelerometer location in front 
of the center of gravity, there is a need for small baffling. The further the accelerometer 
is shifted toward the nose, the less damping is required. The danger zone decreases also 
to one between the center of gravity and center of instantaneous rotation (Figures 10s 
and lot) . 


c. Ideal Accelerometer : For ideal accelerometer control the stability of the 
space vehicle is insensitive to the location of the accelerometer (Figure 9g) . This, how- 
ever, as in all previous cases is valid only for the particular accelerometer gain value 
A = 1. 0. For other gain values A>1, accelerometer location in front of the center of * 
gravity results in instability for the tanks in front of the center of instantaneous rotation, 
while for a location behind the center of gravity all tanks behind the center of instan- 
taneous rotation exhibit instabilities, which decrease with a further shift of the location 
of the accelerometer toward the aft of the vehicle (Figure 9h) . A change of accelerom- 
eter location toward the nose of the vehicle results in less baffling requirements in the 
tanks in front of the center of instantaneous rotation. 

Enhanced stability can always be obtained by changing the previous parameters 
and by proper choice of control values. The sloshing mass can be reduced by special 
selection of tanks for the propellant. Tanks with a large fineness ratio (fluid height over 
diameter ratio) exhibit small sloshing masses. A circular cylindrical tank compared 
with the same tank containing cross walls, thus making four cylindrical quarter tanks, 
exhibits a ratio of the slosh masses of about 1:3. 3. This means that, by introduction of 
quarter tanks, the sloshing mass can be reduced to one third. Besides this the slosh 
frequency increases, thus shifting the natural liquid frequency into a better region (above 
the control frequency: oj g > cu c ) . For four cluster tanks with circular cross sections 
compared with four quarter tanks, the sloshing ratio is still 1:1.6. This means that, by 
the introduction of quarter tanks instead of four cluster tanks of equal volume, the slosh- 
ing mass can be reduced to six tenths. It should be noted, however, that additional 
sloshing masses appear in the quarter tank arrangement, which have only a slightly dif- 
ferent frequency, but can no longer be neglected. These frequencies are about 20% apart. 
The change of the control frequency and control damping is a matter of control system 
design (change of the gain value a Q and a t ) , the selection of center of gravity (which 
reflects on the swivel point location Xg) , the moment of inertia of the vehicle, and the 
longitudinal acceleration of the vehicle. The slosh frequency can be increased by sub- 
division of the tanks in the longitudinal direction. Large accelerometer frequency en- 
hances the stability and the danger zone and thus requires less baffling for stability. 
Proper choice of accelerometer gain values and accelerometer location also plays an 
important role. 
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APPENDIX 


MECHANICAL MODEL 

An equivalent mechanical model which describes the motion of the propellant in 
the tank can be seen in Figure i. This model consists of vibrating mass pointy, springs 
and dampers, a fixed mass with a moment of inertia and a massless disc with a moment 
of inertia. It is so derived that it exerts the same forces and moments as the propellant 
and has the same natural frequency as the liquid. The complete analysis of pre- 
senting the fluid motion as mechanical model either with mass-spring damper systems 
or pendulum systems is given in Reference 3 for cylindrical tanks with circular or 
annular cross sections. Part of the results will be summarized here. 
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B. For a Tank with Annular Cross Section 
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FIGURE 1. MECHANICAL MODEL FOR PROPELLANT SLOSHING 
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JRAL FREQUENCY f FOR LIQUID IN A CIRCULAR 
NDRICAL TANK. n 
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:y f for a liquid in a cylindrical quarter tank 
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FIGURE 4. SLOSHING 




FOR CYLINDRICAL TANK WITH 
ANNULAR CROSS SECTION 
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FOR FOUR CYLINDRICAL QUARTER TANKS 
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FIGURE 6. SI 
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FIGURE 7a. STABILITY BOUNDARY FOR RIGID ACCELEROMETER CONTROLLED 

SPACE VEHICLE (SLOSHING SUPPRESSED) VARYING ACCELEROMETER 
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FIGURE 7b. STABILITY BOUNDARY FOR RIGID ACCELEROMETER CONTROLLED 

SPACE VEHICLE (SLOSHING SUPPRESSED) VARYING ACCE LEROMETER 
DAMPING FOR = 55 rad/sec 




FIGURE 7c. STABILITY BOUNDARY FOR RIGID ACCELEROMETER CONTROLLED 

SPACE VEHICLE (SLOSHING SUPPRESSED) VARYING ACCELEROMTER 
DAMPING FOR co a = 12 rad/sec 
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FIGURE 7d. STABILITY BOUNDARY FOR RIGID ACCELEROMETER CONTROLLED SPACE 
VEHICLE (SLOSHING SUPPRESSED) VARYING CONTROL FREQUENCY 
FOR oo., - 55 rad/sec 




waries \\ INSTABLE 



FIGURE 7e. STABILITY BOUNDARY FOR RIGID ACCELEROMETER CONTROLLED 
SPACE VEHICLE (SLOSHING SUPPRESSED) VARYING CONTROL 
FREQUENCY FOR u - 12 rad/sec 







FIGURE 7f. STABILITY BOUNDARY FOR RIGID ACCELEROMETER CONTROLLED 
SPACE VEHICLE (SLOSHING SUPPRESSED) VARYING CONTROL 
DAMPING FOR w =55 rad/sec 
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FIGURE 7g. STABILITY BOUNDARY FOR RIGID ACCELEROMETER CONTROLLED 
SPACE VEHICLE (SLOSHING SUPPRESSED) VARYING CONTROL 
DAMPING FOR co a = 12 rad/sec 


(NO ADDITIONAL CONTROL) 
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FIGURE 8a. STABILITY BOUNDARY FOR RIGID SPACECRAFT. VARYING 
SLOSH MASS RATIO 






(NO ADDITIONAL CONTROL) 



POSITION GYRO 






(NO ADDITIONAL CONTROL) 
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FIGURE 8c. STABILITY BOUNDARY FOR RIGID SPACECRAFT VARYING 
CONTROL DAMPING 




(NO ADDITIONAL CONTROL) ~i 08 



(NO ADDITIONAL CONTROL) 
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FIGURE 8e. STABILITY BOUNDARY FOR RIGID SPACECRAFT. VARYING 





-1600 -1400 - 1200 -1000 -800 -600 
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FIGURE 9a. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH IDEAL 
ACCELEROMETER CONTROL. VARYING SLOSH MASS RATIO 


-1600 -1400 -1200 -1000 -800 -600 
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FIGURE 9c. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH IDEAL 
ACCELEROMETER CONTROL. VARYING CONTROL DAMPING 









R RIGID SPACECRAFT WITH IDEAL 
3L. VARYING SLOSH FREQUENCY 
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FIGURE 9e. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH IDEAL 

ACCELEROMETER CONTROL. VARYING ATTITUDE GAIN VALUE 







FIGURE 91'. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH IDEAL 

ACCELEROMETER CONTROL. VARYING ACCELEROMETER GAIN 




FIGURE 9g. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH IDEAL 

ACCELEROMETER CONTROL. VARYING ACCELEROMETER LOCATION 
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FIGURE 9h. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH IDEAL 

ACCELEROMETER CONTROL. VARYING ACCELEROMETER LOCATION 
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FIGURE 10a. STABILITY BOUNDARY FOP RIGID SPACECRAFT WITH A REAL ACCELEROMETER 
CONTROL. VARYING SLOSHING MASS RATIO WITH ACCELEROMETER FREQUENCY 
OF 55 RAD/SEC AND GAIN VALUE A = 1. 0 






FIGURE 10b. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH A REAL ACCELEROMETER 
CONTROL. VARYING SLOSHING MASS RATIO WITH ACCELEROMETER FREQUENCY 
OF 55 RAD/SEC AND GAIN VALUE X = 1. 5 
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FIGURE lOd. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH A REAL ACCELEROMETER 
CONTROL. VARYING CONTROL FREQUENCY WITH ACCELEROMETER FREQUENCY 
OF 55 RAD/SEC AND GAIN VALUE A = 1. 0 






68 


FIGURE iOf. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH A REAL ACCELEROMETER 
CONTROL. VARYING CONTROL DAMPING WITH ACCELEROMETER FREQUENCY 




FIGURE lOg. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH A REAL ACCELEROMETER 
CONTROL. VARYING CONTROL DAMPING WITH ACCELEROMETER FREQUENCY 
OF 12 RAD/SEC AND GAIN VALUE A = 1 
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FIGURE 10 h. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH REAL ACCELEROMETER 
CONTROL. VARYING SLOSH FREQUENCY WITH ACCELEROMETER FREQUENCY 
OF 55 RAD/SEC AND GAIN VALUE A = i 
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CONTROL. VARYING SLOSH FREQUENCY WITH ACCELEROMETER FREQUENCY 
OF 12 RAD/SEC AND GAIN VALUE X = 1 






FIGURE lOj. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH REAL ACCELEROMETER 

CONTROL. VARYING ATTITUDE GAIN VALUE WITH ACCELEROMETER FREQUENCY 




FIGURE 10k. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH REAL ACCELEROMETER 

CONTROL. VARYING ATTITUDE GAIN VALUE WITH ACCELEROMETER FREQUENCY 
OF 12 RAD/SEC AND GAIN VALUE A =1 




FIGURE 101. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH REAL ACCELEROMETER 

CONTROL. VARYING ACCELEROMETER GAIN WITH ACCELEROMETER FREQUENCY 
OF 55 RAD/SEC 
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VARYING ACCELEROMETER FREQUENCY 
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CONTROL. VARYING CONTROL DAMPING WITH ACCELEROMETER FREQUENCY 
OF 55 RAD/SEC AND GAIN VALUE A = 1 



« 
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FIGURE iOp. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH REAL ACCELEROMETER 
CONTROL. VARYING CONTROL DAMPING WITH ACCELEROMETER FREQUENCY 





OF 12 RAD/SEC AND GAIN VALUE 







FIGURE lOr. STABILITY BOUNDARY FOR RIGID SPACECRAFT WITH REAL ACCELEROMETER 
CONTROL. VARYING TANK RADIUS WITH ACCELEROMETER FREQUENCY 
OF 55 RAn/SFr 
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NG TANK RADIUS WITH ACCELEROMETER FREQUENCY 









